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Abstract 



An analytical proof of the existence of negative modes in the odd-parity pertur- 
bation sector is given for all known non-Abelian Einstein- Yang-Mills black holes. The 
significance of the normalizability condition in the functional stability analysis is empha- 
sized. The role of the odd-parity negative modes in the sphaleron interpretation of the 
Bartnik-McKinnon solutions is discussed. 
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1 Introduction 



Soon after the discovery of regular particle-like [|l| and black hole solutions in the 
SU{2) Einstein- Yang-Mills (EYM) theory, their instability was demonstrated Nu- 
merical analysis has revealed the existence of n (the number of nodes of the YM function) 
negative modes in the even-parity spherical perturbation sector; moreover, the investiga- 
tion of the non-linear dynamics of perturbations has clearly shown the instability (an 
interesting global analysis can also be found in ||^). Some time later, a generalization 
of the SU{2) black hole solutions for the non- vanishing charge was found in the EYM 
theory for a larger gauge group however, the issue of stability for these solutions has 
been open so far. 

It is worth noting that, for the regular Bartnik-McKinnon (BK) solutions, there 
exists at least one negative mode in the odd-parity sector too. This mode has the same 
nature as negative mode of the electroweak sphaleron solution 0, 0, indicating the 
presence of a potential barrier separating topologically distinct YM vacua. This plays 
the crucial role in the proposed sphaleron interpretation for the BK solutions [jlO|, and 
the related physical issues [|ll|] . The existence of the odd-parity negative mode has been 
shown analytically using variational techniques |10[- The idea consists in the construction 



of an energy-reducing function sequence on the constraint surface in configuration space 
of the theory. The advantage of this technique is that it does not require the detailed 
knowledge of the background equilibrium solution. Recently this approach has been 
successfully applied in the analysis of stability of regular solutions of the EYM-Higgs 
theory |jl2[, and regular solutions of the EYM theory for an arbitrary gauge group |jl3 



The latter result shows, in particular, the instability of the non-trivial SU{3) EYM 



regular solutions found recently by Kiinzle [|T^ 



For EYM black holes, the variational approach can be applied as well. For the 
SU{2) solutions, a sequence of trial functions reducing the total ADM energy has been 
constructed in ||15|. More detailed analysis reveals, however, that the decrease of the 
energy is not sufficient to demonstrate instability. Another important condition is nec- 
essary: the energy reducing perturbations must be normalizable. Unfortunately, for 



the trial functions used in our earlier paper |T^, the normalizability condition was not 
satisfied because of their inappropriate behaviour at the event horizon (for the regular 
solutions such a problem does not arise). Surprisingly enough, we find that a similar 
situation can be often met in other papers, where the argument based on the existence 
of the energy reducing fiuctuations is used — little attention is usually paid to the nor- 
malizability condition. 

The purpose of this paper is two-fold: first, we want to elucidate the nature 
and to emphasize the importance of the normalizability condition for the general field- 
theoretical stability analysis. Secondly, we investigate the structure of the odd-parity 
spherical perturbation sector for the static EYM black hole and sphaleron solutions. We 
shall give a proof of the existence of negative modes in this sector for all known essen- 
tially non-Abelian EYM black holes, and our proof is valid in the regular well. 
This justifies the main statement made in [|15| concerning the SU{2) EYM black holes. 
In addition, we analyse stability of magnetic U{1) black holes within the context of EYM 
theory, and discuss the importance of the odd-parity negative modes for understading 
the physical nature of the localized finite energy EYM solutions. 



1 



2 The functional criteria of instability 



Consider a theory of a classical field, or several fields, denoted commonly by 0. Assume 
that the theory admits a localized static equilibrium solution 0s (x) possessing finite 
energy i?[0s(x)]. Suppose one has found a sequence of static field configurations, 0a (x), 
such that 0A=o(x) = 0s (x), and the energy has a maximum for A = 0: 

E[0A^O(X)] < i?[0s(x)]. (1) 

The mechanical picture for this situation is that of a particle sitting at the top of potential 
hill whose profile is given by U{X) = E[(j)\{x.)]. The question arises whether (|1|) is 
sufficient in order to reveal an instability of the solution 0s (x). Strictly speaking, the 
answer is negative. Namely, the equilibrium state of the "particle" will be unstable only 
if its effective mass is finite. Otherwise, the particle can not depart from the top due to 
infinite inertia. 

Let us reformulate this as follows. If one allows for a time evolution along the 
sequence 0a (x), replacing the parameter A by a function A(t), and if one finds that the 
kinetic energy, i.e. the part of the total energy proportional to A^ , is finite, then (|1[) 
implies indeed an instability. Otherwise, 0a (x) is not a physically acceptable sequence, 
and nothing can be inferred from Eq.(|T]). 

To put this into more rigorous form, consider small fluctuations 5(j){t, x) around 
the static equilibrium solution. Linearizing the field equations and specifying the time 
dependence as (50(t, x) = exp(— zu;t)\E'(x) one can usually represent the perturbation 
equations as []T2[ 

= tu^M^. (2) 

Here, the two operators H and M depending on 0s (x) are assumed to be independent 
on u. H is usually self-adjoint with respect to a properly defined scalar product (^^l^), 
and the kinetic energy matrix M is positive definite, (\I/|M|\E') > 0. 

Clearly, the equilibrium static solution will be unstable if a normalizable solution 
to (D with < can be found. However, this may require numerical calculations. 
The simpler (though less informative) way to reveal an instability is to make use of the 
minimum principle for the following functional defined through (0): 

with \1/ being a trial function. The lowest eigenvalue is known to correspond to the 
lower bound for this functional. From here it follows that if a \1/ can be found such 
that a;^(\I^) < 0, then the operator H is not positive definite, and negative eigenvalues 
therefore do exist. 

Obviously, a;^(^) < imphes that 

<0, (4) 

which is equivalent to (|I|). This condition is frequently used in order to demonstrate in- 
stability in field-theoretical systems (see for instance |^, 0). However, for such systems 
the operator M is generally unbounded. Often, the trial functions from the domain of 
H, while ensuring (4), lead to divergence of the expectation value (\E'|M|\1'). In this case 
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d^) is insufficient to establish instability. We therefore arrive at the second important 
condition 

(^|M|^) < oo. (5) 

Physically this can be seen as the condition assuring the finiteness of the "kinetic energy" 
associated with time-dependent perturbations. 

We deliberately spend so much time on these subtleties since in the existing litera- 
ture they seem often to be overlooked. Sometimes the trial functions used automatically 
fulfill condition but not always. One may think that condition (^) is somehow implied 
by (H), in the sense that, if a function subject to (^) is known, then one can in principle 
find another function satisfying both (H) and (|]). There are, however, examples when 
condition (|) holds even for perfectly stable solutions — but not both (|D and (^. Con- 
sider, for instance, the Bogomolny-Prasad-Sommerfield monopole solutions |T^. They 
form a family 0bps (x,f) whose members are distinguished by the asymptotic value of 
the Higgs field, v, and satisfy the same system of equations; the corresponding mass is 
proportional to v. Given a member of the family specified by some Vq, one can construct 
the trial sequence 0a(x) = (pBPsi'^i'^o ~ -^^) which reduces the energy. However, the 
norm (^) is divergent due to the change of the asymptotic values of the fields. 

Finally, we sketch an idea of the rigorous justification of the above functional criteria 
for instability. Assume that the matrix M in @) is non-degenerate and symmetric. 
Perform a linear transformation ^ = such that O^MO = 1 and the new Hamiltonian 
H = O^HO. Omitting tilde, one can represent (|]) as an eigenvalue problem, = a;^\l/, 
for the unbounded (but usually bounded from below) operator H acting on the Hilbert 
space 7"^ of \I' with a scalar product (\E'|$) = / c?/i(x), where dfi is an appropriately 
chosen measure. Usually, one can easily specify a dense set T>{H) C as the domain 
of H in such a way that H is symmetric on V{H). Next, assume that the existence of 
a self-adjoint extension Hi is known, V{Hi) D V{H). Then, the spectral decomposition 
for Hi implies the following inequality for the ground state eigenvalue: 

u;',<u;\n where o;2(vl/) = i^^M, v]. e (6) 

Thus, if ci;^(\l/) < for some \1/ G V{Hi), then uJq < 0. Notice that if one chooses 
^ e 'D{H) C T>{Hi) then one can replace Hi by H in the definition of uj'^{^). The 
condition < is then equivalent to the following two conditions. The ffist one is 

given by (P, where the trial function \1/ has to be an element of the already specified set 
T>{H) C Ti. The second, the normalizability condition, becomes now fairly trivial 

(^1^) <oo; (7) 

this simply states that \1/ must belong to the Hilbert space 7i. Note that these arguments 
do not assume H to be essentially self-adjoint, which is often not easy to show. It suffices 
to ensure that a self-adjoint extension for H exists, for which one has simple powerful 
criteria |T^. 

In what follows, we apply this procedure within the context of EYM theory. 
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3 Existence of odd— parity negative modes for EYM 
black holes 



All known essentially non-Abelian EYM black hole solutions 0, can be obtained 
within the context of the SU{2) x f/(l) EYM theory. The corresponding action is 

^ = "16^ / ^v^^'^ ' J (i^^^^'^^''' + i-^A^^-^'O ^d'x, (8) 

where F^^, = d^A^ — d^A^ — z[y4^,y4,^] is the matrix valued gauge field tensor, e is the 
gauge coupling constant, A^ = A^r''/2, and r° (a = 1,2,3) are the Pauli matrices; J-'^i, 
is the U{1) field strength. 

In the spherically symmetric case, the spacetime metric is chosen to be 

ds^ = ll {^{l - 2m/r)a^de - ^ - r\dd^ + sin^ M^^)^ . (9) 



Here le = v47r/p;/e is the only dimensional quantity in the problem {Ipi being the Planck 
length); the functions m and a depend on t and r. The U{1) part of the gauge field is 
chosen to be of the dyon type, eJF = (qe/r'^)dt A dr + qmsm{} d'd A dip, satisfying the 
Maxwell equations for any constant qe and qm- 

The spherically symmetric SU{2) YM field can be parameterized by 

eA = WoLi dt + WiLi dr + {p2 L2-(l-pi) L^} di^ + {{l-pi) L2+P2 Lgjsin^? dy^, (10) 

where Wq, Wi, pi, p2 are functions of t and r, Li = n"r"/2, L2 = d^Li, siwdL^ = d^Li, 
and = (sin cos if, sin sin (f, cos . The gauge transformation 

A^ UAU-^ + iUdU-\ with f/ = exp(zfi(t, r)Li), (11) 

preserves the form of the field ([T0|), altering the functions Wo, Wi, pi, p2 as 

Wo^Wo + n, Wi^Wi + n', p± = p^ ± ip2 exp{±in)p±; (12) 

here dot and prime denote differentiation with respect to t and r, respectively. 

It is convenient to introduce the complex variable f = Pi + ip2 and its covariant 
derivative D^f = {d^ — iVF^)/, as well as the (1 + l)-dimensional field strength W^^, = 
dfj.Wi, — dyW^, (/i, z/ = 0, 1). The full system of the EYM equations then reads 

d,{r^aW^^) - 2a Im {fD'^fY = 0, (13) 
D,aD^f-^i\f\'-l)f = 0, (14) 

= + ^1 Wr + N\DJ\^ + ^(l/r - 1)^ + (15) 

m = 2N Re Dof{DJ)\ (16) 

i^--y=l{j^2\Dofn\Dj\'), (17) 

where = 1 — 2m/r, q^ = q^ + q^, and asterisk denotes complex conjugation. 
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Both for g = (the SU{2) case) [Q, and <? 7^ 0, these equations are known to 
possess static non-Abehan black hole solutions labeled by a pair (n, ru), where n is an 
integer, and th > |g| is the black hole radius. For these solutions, Wq = Wi = P2 = 0; the 
function / = Re f = pi has n nodes in the domain rn < r < 00 and tends asymptotically 
to ±1 always remaining in the stripe —1 < /(r) < 1. The metric functions N and a 
increase monotonically from N{rH) = to A^(oo) = 1 and from < u(rH) < 1 to 
cr(oo) = 1, respectively. Asymptotically, the metric is of the Reissner-Nordstrom (RN) 
form with charge q. 

It is worth noting that SU{2) x U (1) group may also subgroup of a larger 

gauge group 0. The basic EYM equations remain of course the same in this case, and 
the only difference may appear in the quantization condition for the Abelian magnetic 
charge g^. 

Consider small perturbations of a given black hole solution 



m 



m + 6m, a^a + 5a, f^f + 6f, Wo = Wi = ^ 5Wo, 6Wi. 



The linearization of the field equation reveals that the even-parity perturbations, 
{6m, 6a, 6pi), and the odd-parity perturbations, {6W0, 6Wi, 6P2), decouple, and therefore 
are independent — this is because the background solutions are invariant under parity. 
Notice that the infinitesimal gauge transformation (|ll|) does not alter the even-parity 
perturbations, while the odd-parity ones change as 



6W0 ^ 6W0 + tl, 6Wi ^ 6Wi + 9!, 6p2 6p2 + 



(19) 



For the even-parity modes, the perturbation equations for {6m, 6a, 6pi) can be 
reduced to a single Schrodinger-type equation. In the g = case, this equation was 
analyzed numerically by Straumann and Zhou who found exactly n bound states for 
any background (n, th) black hole solution. Here we will concentrate on the odd-parity 
sector, and our results do not depend on q. 

For the odd-parity fluctuations, the metric remains unperturbed [12|, thus the 
perturbation equations can easily be obtained via expanding the Yang-Mills equa- 
tions ([T3|) , (0) with respect to 6W0, 6Wi, and 6p2 alone. We use the gauge free- 
dom ( p]P| ) to ensure the condition 6W0 = 0, and specify the time dependence as 
6Wi = ^j2/{r'^N)a{r) exp{—iujt), 6p2 = C,{r) exp{—iujt). The resulting equations can 
be represented as the eigenvalue problem 



where ^E' 




H = aN 



and the Hamiltonian is given by 

/ 2ap/r^ apN/r\f~ifp) \ 

\ (/' + ipf)crpN/r^ paNp + a{p - l)/r'^ y 



(20) 



(21) 



The quantities /, a, N refer to the background static black hole solution, and p = —id/ dr. 
There exists also an additional equation due to the Gauss constraint (Eq. (|1^) with z/ = 0) 



2r2 



a 



(22) 
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One can see however that, as long as 7^ 0, this equation is a differential consequence of 
(HO); (HI)- observe therefore that any solutions to the eigenvalue problem (|20|) (^Tf ) 
with a; 7^ automatically satisfy the Gauss constraint, that is, they have correct initial 
values. 

We introduce the tortoise coordinate r* G R, such that dr^ = dr/{aN), and define 
the inner product as 

/oo roo 
^{^2 = / (al«2 + m)dn. (23) 
-00 J —00 

The Hilbert space can then be chosen to be 7i = L^(R, cir*) © L^(R, dr^). Consider the 
dense set in Ti. consisting of twice continuously differentiable functions with compact 
support. H maps this set into Ti. and it is symmetric on this set with respect to the 
inner product (|23D . We can therefore specify this set as the domain of H, T>{H) = 
Cg(R) X Cq(R) C H 113, (in addition, H can be shown to be bounded from below on 

Notice that H, being real, commutes with complex conjugation. A theorem by Von 
Neumann then ensures the existence of self-adjoint extensions for H. Thus, according 
to the general procedure outlined above, to demonstrate instability it is sufficient to find 
a function G Vi^H) satisfying the inequality (p. 

Consider the set of real functions {hk{r^)} C Cq(R), where hk{r^) = h{r^/k), k > 1, 
and /i(r*) is defined as follows: /i(r*) = h{—r^), h = 1 when r^. G [0,a], —D < h[ < 
for G [a, a + 1], and /i = for r^, > a + 1. Here a, D are positive constants, prime and 
asterisk denote differentiation with respect to r^,, h'^ = dh/dr^, = aNh'. 

Using {hk}, construct the set {"^k} C. 'D{H) specified by 



Vm^a C' ^ ' hk = aohk, ^k = {fi)hk = ^ohk, (24) 



(notice that the background solutions {f, f',N,a) are at least twice differentiable |T^). 
Integration by parts allows us to represent the expectation value {'^klHl'^k) in the form 




{H)= I {{ (6): - faNJ-^ak 1 + 2aNf'J^ak^k + ^{f - 1)4' } dr 



(25) 

Substituting (|2^) into (^5]) and taking into account the background Yang-Mills equation 



= ^f{f - 1) (26) 



we obtain 



/oo (' /2 

dn ^oNcT^ 
-00 

r dneJ{hk)? + 5^{l 

J —00 I T 




nii-hi)}^i+m. (27) 



In this equation, the /c-indepentent first term I is finite and manifestly negative (recall 
that I/I < 1, in addition, ~ — as r^, —00, fl ~ as 00, so the 

integral exists). Taking into account the properties of the smoothening functions hk 
specified above, one can see that the integral I{k) converges uniformly and tends to zero 
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as k increases. We therefore conclude that, for large k, the contribution of I{k) into 
(|27| ) is negligible, thus the whole expression turns out to be negative. This shows that 
with sufficiently large k are such functions from T>{H) that {'^klHl'^k) < 0, which 
completes the proof of the existence of the odd-parity negative modes. Physically, these 
correspond to such admissible time-symmetric initial data which give rise to a growing 
instability [p!9[ . 

Our analysis is valid for any th > \q\- For the uncharged SU{2) solutions, the limit 
r// — > is known to relate to the regular case. Then N{r) > for any r > 0, so r* runs 
over semi-axis, and the Hilbert space is therefore H = L^(R+,(ir*) © L^(R+,(ir*). One 
can see that the analysis presented above is equally valid in this case. 

Finally, for the sake of completeness, we analyse stability of the Abelian magnetic 
f/(l) black holes. For these solutions, the dynamics of perturbations is still governed by 
Eqs. (|20|) , (|2l|) , where the parameters of background solutions are Wo = Wi = / = 0, 
a = I, N = I - 2M/r + = (r - r+)(r - rJ)/r'^, M being the ADM mass. We 
specify T>{H) as before and choose the set {^k} C V^H) parameterized by ak = 0, 
^fc = (r — r^)/r'^hk. Repeating the above procedure, one arrives at 



oo 



oo 



r2^V 420r1 



where I{k) with growing k. This shows that Abelian EYM black holes are unstable 
with respect to non-Abelian fluctuations (see also ||2( 



4 Discussion 

Our analysis reveals the existence of at least one odd-parity negative mode for all 
known non-Abelian EYM black hole solutions, indicating in particular that all of them 
are unstable. For the uncharged SU{2) solutions, taking into account the results of the 
analysis by Straumann and Zhou 0, we therefore conclude that each (n, r^^) black hole 
has at least n + 1 unstable modes - in the even-parity sector and at least one odd- 
parity negative mode. Our analysis is equally valid in the regular case, where such an 
odd-parity mode has precisely the same meaning as negative mode of the electroweak 
sphaleron solution 0, In this sense, this mode is fairly interesting and, from the 
physical point of view, quite typical. In the regular case, it had been precisely the 
existence of this mode which allowed us to suggest a sphaleron interpretation for the BK 
solutions [0. On the other hand, the existence of the even-parity sphaleron negative 



modes is a rather peculiar phenomenon which is present in the EYM theory |rT| . 



For black holes, as has already been mentioned |]T5|, the sphaleron interpretation 
is not so transparent. By definition, a sphaleron is a static solution "sitting" at the top 
of a potential barrier separating topologically distinct YM vacua, whereas in the black 
hole case there are no pure vacuum states because of the finite temperature associated 
with the event horizon (except for the extreme case [^). The existence of the odd- 
parity negative modes in the black hole case shows, however, that the structure of the 
energy surface in the vicinity of the EYM black hole solutions in function space is similar 
to that for the regular BK objects. This suggests one to think of them as "black holes 
inside EYM sphalerons" . An interesting open issue is the exact number of the odd-parity 
modes both for regular and black hole EYM solutions. Investigation of this requires a 
further numerical work. 
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